Influence of quantum fluctuations on zero-temperature phase transitions between 
coflinear and noncoflinear states in frustrated spin systems 
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We study a square- lattice spin-half Heisenberg model where frustration is introduced by competing 
nearest-neighbor bonds of different signs. We discuss the influence of quantum fluctuations on the 
nature of the zero-temperature phase transitions from phases with coUinear magnetic order at small 
frustration to phases with noncoUinear spiral order at large frustration. We use the coupled cluster 
method (CCM) for high orders of approximation (up to LSUB6) and the exact diagonalization of 
finite systems (up to 32 sites) to calculate ground-state properties. The role of quantum fluctuations 
is examined by comparing the ferromagnetic-spiral and the antiferromagnetic-spiral transition within 
the same model. We find clear evidence that quantum fiuctuations prefer coflinear order and that 
they may favour a first order transition instead of a second order transition in case of no quantum 
fluctuations. 



Introduction. While quantum fluctuations do not influ- 
ence the critical properties of phase transitions at T > 
they play an important role at T = and can yield to 
quantum phase transitions, which have attracted a lot 
of attention in recent times (see, e.g., Ref. Quantum 
fluctuations arise due to Heisenberg's uncertainty prin- 
ciple and play a similar role as thermal fluctuations (for 
T > 0) in classical transitions. The spin-half Heisen- 
berg model is a basic model which shows strong quan- 
tum fluctuations in the antiferromagnetic case. While 
the ground-state of the pure Heisenberg antiferromag- 
net (HAFM) on the square lattice shows Neel long- 
range orderi (LRO), a competition of bonds can increase 
quantum fluctuations and may result in rotationally in- 
variant paramagnetic states, suppressing the (coflinear) 
Neel order. This is demonstrated by recent experi- 
ments on (quasi-)two-dimensional Heisenberg systems, 
like CaV409 (see, e.g., Refs. |,|) or SrCu2(B03)2 (see, 
e.g., Refs. 



Besides local singlet formation magnetic frustration is 
an important mechanism to drive zero-temperature tran- 
sitions. In the classical Heisenberg model strong frustra- 
tion often leads to noncoflinear (e.g., spiral) spin states 
which may or may not have counterparts in the quantum 
case. It is generally argued that quantum fluctuations 
prefer a coflinear ordering. A typical example is the frus- 
trated spin-half Ji — J2 model on the square lattice (see, 
e.g., Refs. ^]j8|j9|, p^ , |rT| ) . Here the classical version of the 
Ji — J2 model has a continuously degenerate ground state 
for J2 > but quantum fluctuations can remove 

this degeneracy yielding to a coflinear state ( "order from 



disorder" phenomenon, see, e.g., Refs. |lj,g_3|). More- 
over quantum fluctuations can shift the critical point of 
a coUinear-noncoUinear transition so that the coflinear 
state can survive into a region where classically it is al- 
ready unstable.000 

In this paper we extend our previous work,lll where 
we have studied the transition from a coflinear Neel or- 
der to noncoflinear spiral order in a frustrated spin-half 
HAFM and consider now the transition form a coflinear 
ferromagnetic order to a noncoflinear spiral order within 
the same model. While in the classical version of the 
model both situations can be mapped onto each other the 
quantum model behaves basically different in both cases. 
That is because of the different nature of the coflinear 
state: While the quantum Neel state on two-dimensional 
lattices exhibits strong quantum fluctuations (the sublat- 
tice magnetization of the HAFM on the square lattice is 
only about 60% of its classical value) the ferromagnetic 
state is the same for the quantum and the classical model 
and there are no quantum fluctuations in this state. 

We use the coupled cluster methodEBlii (CCM) and 
exact diagonalization of finite systems to calculate the 
ground state. The CCM is a very powerful method and, 
particularly high-order implementations of this method 
can be used to obtain a consistent description of various 
aspects of quantum spin systems (for an overview see, for 
example, Refs. ^9|j2^,|l|,||j2|||) . We note that another 
important method for spin systems, the quantum Monte 
Carlo method, cannot be used for frustrated spin systems 
since it suffers from the minus sign problem. 

The model. We consider a spin-half Heisenberg model 
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on a square lattice with two kinds of nearest neighbour 
bonds J and J', as shown in Fig. |l|, 

i/ = J^S,.S,+ J'^S.-S,. (1) 

The sums over < ij >i, and < ij >2 represent sums 
over the nearest-neighbour bonds shown in Fig. |^ by 
dashed and sohd hues, respectively. Each square-lattice 
plaquette consists of three J bonds and one J' bond. A 
model with such a zigza g pattern has been treated with 
various methods. 000 

In this paper we consider only the cases in which J 
and J' have different signs (i.e., one bond is ferromag- 
netic while the other is antiferromagnetic) so that the 
plaquettes are frustrated. The case with antiferromag- 
netic J bonds (i.e., J > and hence J' < 0) has been 
studied previously using linear spin wave theory,§l exact 
diagonalization and coupled cluster method0 We there- 
fore focus in this paper our attention mainly on the fer- 
romagnetic case (i.e., J < and J' > 0) but compare 
the obtained results with those of the antiferromagnetic 
case. 

Classical ground state. We consider the ground state 
of the classical version of model i.e., the spins Si 
are assumed to be classical vectors. For |J'| < |J|/3 
(and J and J' having different signs) the ground state 
of (|l]) is coUinear (i.e., ferromagnetic or antiferromag- 
netic depending on the sign of J). At the critical point 
= — J/3, a second-order transition takes place from 
the collinear state to a noncoUinear state of spiral nature 
(see Fig. |l|), with a characteristic pitch angle $ = ±|$ci| 
given by 




Note that for $ = this is the collinear state. 

The spins S^i and Sb, belonging to the A and B sub- 
lattices respectively, can be expressed in terms of the 
spiral k vectoilli with k = (2<I>, 0) (and see Fig. I). We 
note that this spiral state is incommensurate in the a;- 
direction. We also note that for the classical model the 
antiferromagnetic case can be transformed into the ferro- 
magnetic case by the simultaneous substitution J — > — J, 
J' — J', SigB — SigB. Hence the physics for both 
cases is classically the same. 

Calculation of the quantum ground state. To calculate 
the quantum ground state of the Hamiltonian (^ we use 
the CCM. Details concerning the treatment of the model 



d) with the CCM are given in Ref. |T|. We use the CCM 
for high orders of approximation up to LSUB6 (using 
1638 fundamental configurations). 

We further exactly diagonalize finite lattices of rectan- 
gular shape (La: x = 4 x 4, 6x4, 8x4) using periodic 
boundary conditions. The longer side of the rectangle 
corresponds to the direction of the J' bonds and so we 
can diminish the influence of the boundary conditions by 
an increase of . 

The collinear-noncollinear transition. While classi- 
cally we have always a second-order phase transition from 
collinear order to noncoUinear order at = — J/3, we 
obtain for the quantum case a different behaviour for the 
ferromagnetic and the antiferromagnetic case. 

Using the CCM we find for the antiferromagnetic case 
( J = +1) indications for a shift of this critical point to a 
value J^ « —1.35 (see Fig. ||). On the other hand for the 
ferromagnetic case (J = —1) we do not find such a shift 
(see Fig. ^). The exact diagonalization (ED) data of the 
structure factor S'(k) (see Figs. ^ and ||) agree to these 
findings. For J = -1-1 (see Fig. ||) the collinear Neel order 
[k = (0, 0)] becomes unstable against the noncoUinear 
spiral order [k — (7r/4, 0)] in the classical model for J' < 
—0.36 but in the quantum case only for J' < —0.95. The 
situation for the ferromagnetic case (J = —1) is again 
different. Here the results of the structure factor (see 
Fig. H) show that the transition from k = (0, 0) (collinear 
ferromagnetic order) to k = (7r/4, 0) (spiral order) takes 
place at nearly the same value of J' w 0.36 for both, the 
classical case and the quantum case. 

Taking the deviation of the on-site magnetic moment 
{Si) from its classical value {Si)c\ = 1/2 as an indication 
for the degree of quantum fluctuations we can compare 
the strength of quantum fluctuation near the collinear- 
noncollinear transitions for both, the antiferromagnetic 
and the ferromagnetic case. As reported in Ref. |l^ for 
J = the quantum fluctuations are particularly strong 
near the antiferromagnetic-spiral transition leading to an 
on-site magnetic moment less then 20% of its classical 
value. On the other hand, it can be seen from Fig. ^ 
that the on-site magnetic moment takes its classical value 
1/2 up to J' « 0.36 for J = — 1 and therefore virtu- 
ally no quantum fluctuations occur at the ferromagnetic- 
spiral transition. Hence the shift of the critical J[. in the 
antiferromagnetic case can clearly be attributed to the 
strong quantum fluctuations. In general our flndings are 
consistent with the statement that quantum fluctuations 
(which we have in the antiferromagnetic case only) pre- 
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fer a coUinear ordering, so that in this case the quantum 
colhnear state can survive into a frustrated region where 
classically the collinear state is already unstable. 

We further note an agreement between the CCM re- 
sults and the ED results beyond the critical J^. By 
examining the structure factors (see Figs. || and ||) we 
find that for the antiferromagnetic (respectively ferro- 
magnetic) case the transitions to a spiral state with a 
greater k vector (i.e., with a greater pitch angle <&) occur 
in the quantum model always at a greater (respectively 
smaller) absolute value of J' then the corresponding clas- 
sical transitions. This agrees with the CCM results of 
the pitch angle (see Fig. ||), where we have <i>qm < *&ci 
(respectively > ^d)- 

The discussion given above corresponds to our finding 
concerning the order of the transition. Clearly in the 
ferromagnetic case (J' — —1) both the classical and the 
quantum model show a second-order phase transition (see 
Figs. I and I). On the other hand it was discussed in 
Ref. ^ that the coUincar-noncoUinear transition in the 
antiferromagnetic case (./ = +1) is probably a first-order 
phase transition for the quantum model (c.f.. Fig. ^) in 
difference to the second-order transition in the classical 
case. 

Formation of local singlets. For sufficient strong anti- 
ferromagnetic J' bonds the system (|^) is characterized by 
a tendency to singlet pairing of the two spins coupled by 
a J' bond and hence the long-range magnetic (collinear 
or noncoUinear) order is destroyed. We obtain clear indi- 
cations of a second-order phase transition to a quantum 
paramagnetic dimerized phase at a certain critical value 
of J' = J'g. For J' > Jg the on-site magnetic moment 
{Si) becomes zero. For the ferromagnetic case (J — —1) 
we find J'g K, \ using the extrapolated CCM-LSUBn re- 
sults (see Fig. ||). We may also consider the inflection 
points of {Si) versus J' for the LSUBn approximations, 
assuming that the true curve will have a negative curva- 
ture up to the critical point. We find the corresponding 



infiection points at J' « 1.2 (n = 2), J' « 0.76 [n = 4) 
and J' « 0.74 (n=6), indicating a critical even slightly 
smaller then w 1. 

The ED data give a similar approximation of J^: For 
J' « 1 finite-size effects in spin-spin correlations disap- 
pear almost completely (see for illustration Fig. ||). This 
indicates that spin-spin correlations are short ranged 
with a length scale smaller than the size L^. 

We note that for the antiferromagnetic case (J = +1) 
the strength of antiferromagnetic J' needed for breaking 
Neel order by formation of local singlets is much larger 
( « 3, see Ref. |l^). The lower critical ~ 1 in the fer- 
romagnetic case is due to frustration which assists local 
singlet formation (c.f., Ref. |27| , p8| ). 

Summary. Using the CCM and the ED technique 
we have studied the influence of quantum fluctuations 
on zero-temperature transitions between collinear or- 
dered and noncoUinear ordered states in a frustrated 
spin-half square-lattice Heisenberg model with two kinds 
of nearest-neighbour exchange bonds. The frustra- 
tion drives a second-order transition between collinear 
(antiferro- or ferromagnetic) and noncoUinear (spiral) 
states in the classical model. For the quantum model 
the CCM provides a consistent description of collinear, 
noncoUinear, and disordered phases, while some other 
standard techniques (e.g., QMC) are not applicable. We 
find a strong infiuence of quantum fluctuations on the 
nature of the coUinear-noncoUinear transition, and quan- 
tum fluctuations (which favour collinear ordering) may 
change the second-order classical transition to a flrst- 
order quantum transition. If quantum fluctuations are 
suppressed in the collinear phase of the quantum model, 
the transition to the spiral phase is similar for the quan- 
tum and for the classical model. 
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FIG. 1: Illustration of the classical spiral state for the 
square-lattice Heisenberg model of Eq. (|l|, with two kinds 
of regularly distributed nearest-neighbour exchange bonds, J 
(dashed lines) and J' (solid lines). The spin orientations at A 
and B lattice sites are defined by the angles 9n = n$ where 
n — 0,1,2,..., and "I> is the characteristic angle of the spiral 
state. The state is shown for "I> = 7r/12 and n = 0, 1, ... ,7 and 
refers to the ferromagnetic case (J < 0) with a J' > |J|/3. 
For the antiferromagnetic case (J > and J' < —J/3) all 
spins on the B sublattice are reversed. 





FIG. 2; Ground-state structure factor S'(k) cx 
X/i jeA e'''^^ ''(Si . Sj) (i.e., the summation is taken over 
one sublattice) for a 8 x 4 lattice (with antiferromagnetic 
J = -1-1) for the quantum and the classical case for various 
spiral vectors k. 
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FIG. 3: Ground-state structure factor ^(k) (see Fig. ^ for a 
8x4 lattice (with ferromagnetic J — —1) for the quantum 
and the classical case for various spiral vectors k. 
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FIG. 4; Pitch angle <1> versus \ J'\ for the quantum and the 
classical case. While "1> is classically the same for the ferro- 
magnetic case (J = — 1, J' > 0) and for the antiferromagnetic 
case (J — +1, J' < 0) the quantum pitch angle is different 
for both cases. The curves left of the classical (dashed) curve 
belong to J = — 1 and those right of it to J = -1-1. Note that 
for the ferromagnetic case (J = —1) for J' > 1 the pitch an- 
gle "1> becomes meaningless, since the spiral order is already 
destroyed in this region (c.f., Fig. n). 
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FIG. 5: On-site magnetic moment versus J' for the ferro- 
magnetic case {J — —1) calculated within the CCM-LSUBn 
approximations and extrapolated^ to n — oo (The extrapola- 
tion is done as described in RefJij). 




FIG. 6: Nearest-neighbour spin-spin correlation of the two 
spins connected via a J' bond versus J' for the ferromagnetic 
case using exact diagonalization (ED) data. 



